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Abstract 

We present a higher order Generahzed (Gravitational) Uncertainty Principle (GUP) in the form 
[X, P] = ih/{l — /3-P^). This form of GUP is consistent with various proposals of quantum gravity 
such as string theory, loop quantum gravity, doubly special relativity, and predicts both a minimal 
length uncertainty and a maximal observable momentum. Here we present two equivalent and 
exact representations where one is formally self-adjoint and naturally perturbative. We show that 
the presence of the maximal momentum results in an upper bound on the energy spectrum of 
the free particle, particle in box, and harmonic oscillator. We then generalize this GUP to D 
dimensions that will be shown it is noncommutative and find the invariant density of states. This 
form of GUP modifies the blackbody radiation spectrum at high frequencies and predicts a finite 
cosmological constant. Although it does not solve the cosmological constant problem, it gives a 
better estimation with respect to the presence of just the minimal length. 
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I. INTRODUCTION 



In recent years, there is a great interest to study the effects of the Generahzed Uncertainty 
Principle (GUP) and the Modified Dispersion Relation (MDR) on various quantum mechan- 
ical systems. Indeed, the ideas of GUP and MDR arise naturally from various candidates 
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of quantunaffravity such as string theory [lHJ|, loop quantum gravity [5[, noncommutative 
spacetime [6|-|8|], and black holes gedanken experiments BQ. These tKeone.,„d,cate that 
the Heisenberg uncertainty principle should be modified to incorporate additional constraints 
in the presence of the gravitational field. 

The existence of a minimal length scale of the order of the Planck length ipi = ~ 
10~^^m is one of the main outcomes of various GUP proposals where G is Newton's grav- 
itational constant. In fact, beyond the Planck energy scale, the effects of gravity are so 
important which would result in discreteness of the very spacetime. Notably, the quantum 
field theory in curved background can be renormalizable by introducing a minimal observ- 
able length as an effective cutoff in the ultraviolet domain. Also, in the string-theoretic 
argument, we can say that the string cannot probe distances smaller than its own length. 

The introduction of this idea has drawn much attention in the literature to study the 
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effects of GUP on small scale and large scale systems [1 lH26| . It is also possible to incor- 
porate the idea of a maximal observable momentum into this scenario. In fact, in doubly 
special relativity (DSR) theories, we consider the Planck energy (Planck Momentum) as an 
additional invariant other than the velocity of light 27H29|. Recently, the construction of a 
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34|. 



perturbative GUP which is consistent with DSR theories is also discussed in Refs. 
It is also shown that a minimum uncertainty in momentum can arise from curvature, as part 
of a study that indicated that curvature and noncommutativity can be seen as dual to each 



other 



35| 



In this paper, we investigate the effects of a new generalized uncertainty principle to all 
orders in the Planck length on some quantum mechanical systems. This form of GUP implies 
the existence of a minimal length uncertainty and a maximal momentum in agreement with 
various theories of quantum gravity. We present two equivalent exact representations of 
this modified algebra where the latter is self-adjoint and does not change the nature of the 
position operator. Moreover, the existence of the maximal momentum is manifest from this 
representation. 



Here, we study the problems of the free particle, particle in a box, and the harmonic 
oscillator in this scenario and obtain their quantum mechanical and semiclassical energy 
spectrum which are bounded from above. We also address the maximal localization states, 
generalization to D dimensions, validity of semiclassical approximation, invariant density of 
states, cosmological constant, and blackbody radiation in this GUP framework. 



II. THE GENERALIZED UNCERTAINTY PRINCIPLE 

First let us consider a generalized uncertainty principle proposed by Kempf, Mangano 
and Mann (KMM) and results in a minimum observable length 

AXAP>^(1 + /3(AP)2 + C), (1) 

where /3 is the GUP parameter and C is a positive constant that depends on the expectation 
values of the momentum operator, i.e., ^ = We also have {3 = Pq/^Mpic)"^ where 

Mpi is the Planck mass and /3o is of the order of the unity. It is straightforward to check 
that the inequality relation ([T]) implies the existence of a minimum observable length as 
(AX)J^^^ = hy/p. In one-dimension, the above uncertainty relation can be obtained from 
the following deformed commutation relation: 

[x,P] = ih{i + pp^). (2) 

As KMM have indicated in their seminal paper, we can write X and P in momentum space 

Q 

representation as [3J 

P0(p)=p0(p), (3) 
X<P{p) = in{l + l3v')dp<i){p), (4) 

where X and P are symmetric operators on the dense domain 6*00 with respect to the 
following scalar product: 

where i+^p^ \p){p\ = 1 ^ind {p\p') = (1 + 6{p — p'). With this definition, the com- 
mutation relation ([2]) is exactly satisfied. 



Based on the field theory on nonanticommutative superspace, Nouicer has suggested the 
following higher order GUP which agrees with ([2]) to the leading order and also predicts a 
minimal length uncertainty 



[X,P] = ihexp {f3P^) . 



(6) 



This algebra can be satisfied from the following representation of the position and momentum 
operators: 



X(f){p) = ihexp [(3p'^) dp(f){p). 



(7) 
(8) 



Now the symmetricity condition of the position operator implies the following modified 
completeness relation and scalar product 



+00 



J — oo 

{p\p) = exp 6{p- p). 
Also, the absolutely smallest uncertainty in position is given by (AX) 



(9) 
(10) 



Nouicer 



To incorporate the idea of the maximal momentum, Ali, Das and Vagenas have proposed 



the following modified commutation relation 



32|-|34| 



a P6. 



^]+a^ {P% + 3P.P,) 



(11) 



where a = ao/Mpic = aoipi/h is the GUP parameter, = ^ PjPjy ^pi is the Planck 

i=i 

mass, and Mpic^ ~ lO^^GeV is the Planck energy. This form of GUP implies both a minimal 



length uncertainty and a maximal momentum uncertainty, namely 

AX > {AXU. 
AP < (AP)„,. 



32| 



Mpic 



"0 



1/a. 



(12) 
(13) 



The commutation relation ffTTj) is approximately satisfies by the the following representation 



Xi 



Xi 



Pi = Pi (^1 — ap + 2a^p^^ 



(14) 
(15) 



where Xi and pi obey the usual commutation relations [a;i,pj] = iTiSij and p is the magnitude 
of p. Now Eq. ( IT2|) implies a ~ a/5- However, this proposal has the following difficulties: 



It is perturbative, i.e., it is only valid for small values of the GUP parameter. 



• Although the minimal length uncertainty can be considered in some sense as the 
minimal length, the maximal momentum uncertainty cannot be interpreted as the 
maximal momentum in the same way which is required in DSR theories. Indeed 
Eq. f|T3l) puts an upper bound on the uncertainty of the momentum measurement, not 
on the value of the observerd momentum. 

• It does not imply noncommutative geometry, because [Xj,Xj] = [see Eq. ([H])]. 

To overcome these problems, consider the following higher order generalized uncertainty 
principle (GUP*) which implies both the minimal length uncertainty and the maximal ob- 
servable momentum 

This commutation relation agrees with KMM's and Noucier's proposals to the leading order 
and contains a singularity at = This fact shows that the momentum of the particle 
cannot exceed \/a which agrees formally with Eq. (fT3|) . As stated before, Eqs. ( !T3|) 

and f|T6l) imply two basically different quantities. However, the presence of an upper bound 
on the momentum properly agrees with DSR theories. As we shall see, the physical observ- 
ables such as energy and momentum are not only nonsingular, but also are bounded from 
above. 

To satisfy the above commutation relation, we can write the position and momentum 
operators in the momentum space representation as 

P0(p)=p0(p), (17) 

= r^^p'^(^)- (18) 

Using the symmetricity condition of the position operator the modified completeness relation 
and scalar product can be written as 

/•+l/v^ 

{m= / dp{i-pp')r{p)Hp), (19) 

^ (20) 
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TABLE I. The minimal length uncertainties and maximal momentums in three GUP frameworks. 





KMM 


Nouicer 


GUP* 








1 

7^ 


P 

^ max 







The uncertainty relation that arises from GUP* is given by 

h/2 



(AX)(AP) > 



1 -/3P2 



>^{l + P{P^) + P^P^) + P^P^) + ...) 
>\(l + P {P') + {P'f + {P'f + ■ ■ ■ 

> ^ (l + /3 [(AP)2 + {Pf] + [{APf + {Pf] ' + [(AP)2 + (P)2] ' + 

> (21) 

- 1-/3[(AP)2 + (P)2]' ^ > 

where we have used the property (P^") > (P^)". In order to find the minimal length 
uncertainty of this deformed algebra, we consider the physical states for which we have 
(P) = and solve the following saturate GUP* for AP 

<^^)(^^' = t^Utt- '''' 

which has a minimum at AP = 1/ ^/3^. So the absolutely smallest uncertainty in position 
is given by 

(A^)™n = ^^V^- (23) 

In Table [H we have compared minimal length uncertainties from various GUP scenarios. 
These results show that (AX)Kmm < (AX)^™ < (AX);,,,. 



III. FUNCTIONAL ANALYSIS OF THE POSITION OPERATOR 

The eigenvalue problem for the position operator in the GUP* framework and in the 
momentum space is given by 

ih 



1 -/3p2 

6 



apV^A(p) = A^a(p). (24) 
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FIG. 1. Plotting (V'aIV'a') over A - A' in units of fi^ for the KMM GUP (blue line) and GUP* 
(red line). 



This equation can be solved to obtain the position eigenvectors 



V-aIp) = cexp ( ^ M - ) ). (25) 



The eigenf unctions are normalizable 



, "^/^ 4cc* 
1 = cc* / dp (1 - = -j= (26) 



Therefore 



a/3a/^ f-i\p /3 2 



V^a(p) = exp \^i-^p^]]. (27) 

Now we calculate the scalar product of the position eigenstates 



2(A-A') V 

Thus, similar to the KMM scenario, the position eigenstates are generally no longer orthog- 
onal. In Fig. [H we have depicted (ipxlipy) for the KMM proposal and GUP*. Although 
this quantity in both models has a same functional form, it is more oscillatory in the KMM 
framework. 
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A. Maximal localization states 



The maximal localization states ) are defined with the properties 



and 



These states also satisfy 



(AX): 



0. 



(29) 



(30) 



(31) 



2(AP)2 

To proceed further we need to express {[X,P]) in terms of AP and (P). However, since 
{[X,P]) also depends on (P^), {P^), etc., and these quantities cannot be calculated before 
specifying \^|J), to first order in the GUP parameter we can use the approximate relation 
{[X, P]) ih{l + /3(AP)^ + f3{P)'^). So, in momentum space, the above equation takes the 
form 



ih d 
1 — f3p'^ dp 

which has the solution 



(X) + ih 



l + l3{APf + (3{Py 
2(AP)2 



p-(P))')^(p)^0, 



(32) 



ip{p) ~A/'exp 



P 3 



1^^^ + 2(APp 

1 + /3(AP)2 + /3(P)2 V-^p^ 



4(AP)2 V" 2" 7 J' ^^^^ 

To find the absolutely maximal localization states we need to choose the critical momentum 
uncertainty AP = 1/a/3^ that gives the minimal length uncertainty and take (P) = 0, i.e.. 



7/'fL(p) -A/'exp 



-(3 P 



P 4 

2^ 



(34) 



where the normalization factor is given by 



1 = MM* I dp (1 - exp (2/3p2 _ /3y) 



1.0123 



(35) 



Note that ip^^ij)) exactly satisfies Eq. (!29|) . However, because of the approximation that 
assumed to find ip^^{p) fl34l) . it approximately obeys relation fl30l) . i.e.. 



AX|^ML^ 



1.0998(AX); 
8 



(36) 



which shows an error less than 10%. Also, because of the fuzziness of space, these maximal 
localization states are not mutually orthogonal. It is worth to mention that, in this frame- 
work, the expectation value of the kinetic energy operator P"^ /2m is finite for both \ipx) and 
Indeed we have 

p2 



and 



ML 



2m 



p2 



2m 



ML 



1 

lOm/3' 



0.7345 
lOm/3 ■ 



(37) 



(38) 



These quantities for the KMM proposal are oo and l/2m/3, respectively. 
To find the quasiposition wave function iplC,), we define 



ML I 



(39) 



where in the limit /3 — > it goes to the ordinary position wave function = (■Cl^)- Now 
the transformation of the wave function in the momentum representation into its counterpart 
quasiposition wave function is 



m = A/" 



dp (1 — exp 



o I 2 4 

P [P - 2 P 



Hp)- (40) 



This relation shows that similar to the ordinary quantum mechanics and the KMM proposal, 
the quasiposition wave function of a momentum eigenstate ij^p{p) = S{p — p) with energy 
E = p^ /2m is still a plane wave but with a modified dispersion relation 

2T,h Ao,d(^) 



y/2^ (1 - |m/3E) 1 - ^m(3E 



(41) 



where Aord(-E') = 2TTh/\/2mE is the wavelength in the absence of GUP. In Fig. |2]we have 
depicted A versus mE in various scenarios. Since Eq. ( HTl) is bounded from below, there 
exists a nonzero minimal wavelength. So the wavelength components smaller than 



Ar 



3 



(42) 



are absent in the Fourier decomposition of the quasiposition wave function of the physical 
states. Therefore, the maximal energy of a momentum eigenstate is 

3 



E„ 



2ml3' 



(43) 



HE) 




mE 



1 2 3 4 5 6 7 



FIG. 2. The wavelength of the quasiposition wave function of a momentum eigenstate in ordinary 
quantum mechanics, KMM's GUP and GUP* for (3 = 0.2. 



Since the transformation (l40l) as the generahzed Fourier transformation is invertible, the 
transformation of a quasiposition wave function into a momentum space wave function is 
given by 



A/" 



-1 /• + 00 



2Trh 



exp 



n I 2 P 4: 



exp 



h 



m- (44) 



Now the scalar product of states in terms of the quasiposition wave functions reads 



2nh 
X exp 



/ / dpdede'(l-/3/)exp 

-'-oo J-oo 



1/v^ J —CO J—OD 



(45) 



IV. HARMONIC OSCILLATOR 

In this section, we apply the developed formalism to the case of a linear harmonic oscil- 
lator. Using the expression for the Hamiltonian 
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and the representation for X and P, we obtain the following form for the stationary state 
Schrodinger equation: 

d'^^ip) 2(3p di){p) 



+ (l-/3/)2(e_^y)^(p) = 0, 



dp"^ 1 — Pp"^ dp 
where —1/ y/P < J5 < 1/ and 

2E 



(47) 



(48) 



A. The quantum mechanical solution 



Using the dimensionless variable u = \^p, Eq. fj47|) can be written as 



dV(u) ^ 2u di){u) 
dv? 1 — du 

where — 1 < -u < 1 and 



+ (1-M2)'(e'-VV)^(M) =0, 



Now by changing the variable to a; = m — (1/3)m^ we have 



^^^ + r/'2^(x)^(x) = eV(x), 
dx^ 



where —2/3 < a; < 2/3 and 



V{x) 



1 - iVS + (-2)1/3 (3s + V9x2 - 4) 



2/3 



22/3 (3x + V9x2 -4)'^^ 
is the effective potential (see Fig. [3]). The boundary condition now is 



0. 



(49) 



(50) 



(51) 



(52) 



(53) 



To solve Eq. (15T|) . we can expand the wave function in terms of the particle in a box 
eigenfunctions. Since the potential term V{x) is an even function of x, to avoid large 
matrices, we use 



COS 



1\ Tlx 



m 



2 L 



and 



C(^) 



/I /rmrx\ 



(54) 



(55) 
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V{x) 

1 h 




FIG. 3. The effective potential V{x) (blue line) and the ground state wave function (red line) for 
r?' = 100. 



basis functions (m = 1, 2, . . .) for even and odd parity solutions, respectively, and write the 
wave function as il){x) = Ylm ^m4'm{x) which vanishes at ±L. Now the boundary condition 
dMD reads L = 2/3. 

The approximate solutions are the eigenvalues and the eigenfunctions of the (A^ x N) 
Hamiltonian matrix Hat in the form 



m — 



1 \ 2 2 
L2 



+ DI 



and 



Hry 



D 



odd 



(56) 



for even and odd states, respectively. Here, 6mn is the kronecker's delta and 

V2 rL 



Dl 



T] 



V{x) cos 



m 



nx 

T 



cos 



n 



TTX 

T 



dx, 



odd 



L 



(57) 
(58) 



// im^X\ /T17XX\ 
V{x) sin y j si'^ v~^/ 

where m and n run from 1 to N . In the usual diagonalization scheme with the particle in 
box basis functions, we need to adjust the domain L with respect to the number of basis 
functions in such way that the total error to be minimized 36|. However, for our case, 
since the boundary condition flS^ has fixed the domain, i.e., L = 2/3, the accuracy of the 
solutions grows as the number of the basis increases. In Table [TTl we have reported the first 
ten energy eigenvalues of the harmonic oscillator in the GUP* framework. Indeed = 30 
basis functions suffices to obtain nearly accurate results for the low lying energy eigenstates. 
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TABLE II. The energy eigenvalues of the harmonic oscillator in the GUP* framework. Here 
£^ = e'Jr]' = e„/?7 = 2Enlhu, N = 30, and r]' = 100. 



1 

3 
5 
7 
9 
11 
13 
15 
17 
19 



I. 00251 
3.02284 
5.06411 
7.12698 
9.21216 

II. 3204 
13.4524 
15.6091 
17.7913 
20.0000 



I. 00509 
3.02559 
5.06704 
7.13011 
9.21550 

II. 3240 
13.4563 
15.6133 
17.7958 
20.0049 



2.6 
9.1 
5.8 
4.4 
3.6 
3.2 
2.9 



X 10-3 
X 10-4 
X 10-4 
X 10-4 
X 10-4 

X 10-4 

X 10-4 



2.7 X 10-4 
2.5 X 10-4 
2.4 X 10-4 



B. The semiclassical solution 



The total energy in terms of ordinary variables is 



2 2 2 

E^-^ + 2. 59 

To find the approximate energy eigenvalues of the above Hamiltonian, we use the Wilson- 
Sommerfeld quantization rule in the form 

^xdp^(^+^h, 71 = 0,1,2,..., (60) 

where we have used ^ d{xp) — — fxdp + ^pdx. This integral can be written as 

/xdp— / (l — (3p^) \J — p^ dp, (61) 



where z — \j2mE. So the semiclassical energy spectrum is given by 
1 - 2m^nuj in + \) 



(63) 
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10 20 30 40 



FIG. 4. Comparing semiclassical harmonic oscillator energy spectrum for KMM GUP (red line) 
and GUP* (blue line). We set 7 = 0.01. 



where 7 = Pmhu = 77'"^. As it is shown in the appendix, the first three terms are similar 
to the energy spectrum of the harmonic oscillator in the KMM framework. In Fig. HI we 
have depicted the energy spectrum in both GUP^^^ and GUP* frameworks. Note that, 
in GUP* scenario, the energy is also bounded from above. Indeed, the maximum possible 
energy for the harmonic oscillator is 

and the number of states {N = n + 1) is finite, namely 

(65) 



1 1 

,27 ~ 2 



where \_x\ denotes the largest integer not greater than x. So, to have at least one state, we 
should have 

7 < 1. (66) 

As Fig.Hlshows, we have 50 states for 7 = 0.01. The first ten semiclassical energy eigenvalues 
are presented in Table HTl js^]- As the table shows the semiclassical results agree well with 
the quantum mechanical energy spectrum. In fact, the relative error is less than 3 x 10~^ 
even for the ground state. 
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C. The classical solution 



In the classical domain, the equations of motion are 

X = {X,H} ^ 



P = {P,H} 



Tn{l- I3P^)' 
mu'^X 



The solutions to these equations are 



= ( 1 - I) arccos - ^P{tW Pl^^ - P\t), 



X{t) 



l-(3P\t) dP{t) 



where 



V2mE. 



To first-order in /3 we have 

P(t) Pma 
X{t) = Xrru 



COS 



1 + 2 ) 



1 + - ) cos 



— - sin ut cos ut 



- sin ujt 
2 



1 + - 1 



(67) 
(68) 

(69) 
(70) 

(71) 

(72) 
(73) 



where X, 

to C(/3) 



\/ IE I rauj'^ . As we have expected these results agree with the KMM proposal 



It is straightforward to show that the infinitesimal phase space volume between equal 
energy contours E and E + d£', and equal time contours t and t + dt can be written as 



dEdt = (1 -/3p2)dXdP. 



(74) 



Now, since by definition the left hand side of this equation is time independent, the right 
hand is also time independent. 



V. A FORMALLY SELF- ADJOINT REPRESENTATION 

Although the set of Eqs. ( |T71) and (|T8l) is an exact representation of the algebra presented 
in Eq. ( !T6|) . it is not self-adjoint and does not preserve the ordinary nature of the position 
operator. To solve these difficulties, we can write P = /(p) and preserve the ordinary nature 
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df 1 

of the position operator, i.e., X = x. Thus, using Eq. ffT6ll we find — = -— 77 which 

dp 1 - /3f^ 

resuhs in 

f{p)-lf3f{p)=p. (75) 

Consequently, the alternative representation in exact and perturbative form can be written 

as 

X = x, (76) 

/ , N 2/3 

(-2/3)1/3 (sp + ^gp2 _ 4/13) 

P = \ M/3 ^ > (77) 

(2/3)2/3 (3j9+ v/V-4//3 



= P+ \Pp' + V + V + ■ ■ • ■ (78) 

Note that this representation is formally self-adjoint, i.e., A = ioi A E {X, P}. Also 
Eq. fl77|) shows that the momentum P has a upper bound 

2 

at p = Now X and P are symmetric operator on the dense domain 5*00 with respect 

to the following scalar product in the momentum space: 

2 

r{p)<p{p)dp. (80) 

We have schematically depicted the behavior of P versus p in Fig [51 

In this representation, to write the Hamiltonian, it is more appropriate to use Eq. (ITHj) 
and express the Hamiltonian perturbatively as 

H = ^ + V{x) + + + P^^^ + 0(/3^), (81) 

2m Sm 18m 9m 

rn 

which agrees with perturbative version of the KMM proposal to 0{f3) |23i]. In the quan- 
tum domain, this Hamiltonian results in the following generalized Schrodinger equation in 
quasiposition representation: 

- - - - - -(^*) - >'(^)^'<^) - -^•(^'■(-) 

where the extra terms are due to the GUP-corrected terms in Eq. ( IHT]) . Notice Eq. ( 1771) 
shows that this representation is naturally perturbative that is apparent from Eq. (l82l) . 
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V7 





FIG. 5. Schematic behavior of P versus p in the second representation for the ordinary quantum 
mechanics (red Une) and GUP* (blue Une). 



Note that for an operator A which is "formally" self-adjoint {A = A'^) such as (1761) 
and (178|) . this does not prove that A is truly self-adjoint because in general the domains 
T>{A) and P(y4"'') may be different. The operator A with dense domain T>{A) is said to be 
self-adjoint if 'D{A) = V{A^) and A = A"'". For instance, in this representation, X (1761) is 
merely symmetric but not self-adjoint. To see this note that in this representation and in 
the momentum space the wave function 0(p) have to vanish at the end of the p interval 
{—2/3^/]3 < p < 2/3a/^), because the tangent function diverges there. So X is a derivative 
operator ihd/dp on an interval with Dirichlet boundary conditions. But this means that X 
cannot be self-adjoint because all candidates for the eigenfunctions of X, (the plane waves, 
which are even normalizable) are not in the domain of X because they do not obey Dirichlet 
boundary conditions. Calculating the domain of the adjoint of X shows that it is larger 
than that of X, so X is indeed not self-adjoint, i.e.. 



^3V? / d 

dpip*{p) {i^-K- 

37? V ^ 



"ay? 



+ihip*ip)4>ip) 



ihip* {p)(p{p) 



2 



(83) 



3v^ 



Now since (f){p) vanishes at p 



ip* (p) can take any arbitrary value at the boundaries. 



The above equation implies that X is symmetric, but it is not a self-adjoint operator. 
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Although its adjoint = ihd/dp has the same formal expression, it acts on a different 
space of functions, namely 

^ . ^) ; . (^) - . (5^) ^ 0}, (84) 



= 1^/','?/'' £ (^3y^' 3v^} ' other restriction on tp^. 



(85) 



On the other hand, since there is no Dirichlet boundary conditions on the wave functions in 
the position space (—00 < a; < 00), P is still self-adjoint. 

To check the self-adjointness property of X, we can also use the the von Neumann's 
theorem 
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40|. First we need to find the wave functions that satisfy the eigenvalue 
equation 

X^Mp) = ^fidpMp) = ±^A0±(a;), (86) 

that reads 

Mp)=C±e^^P. (87) 

Since both (f)±{p) belong to (3^' 3^)' deficiency indices are (1, 1). Therefore, the 
position operator is not self-adjoint but has a one-parameter family of self-adjoint extensions 
which is in agreement with the previous result. 

A. Free particle 

In ordinary quantum mechanics, the free particle wave function Up{x) is defined as the 
eigenfunction of the momentum operator, namely 

Pup{x) = p Up{x) , (88) 

where p is the eigenvalue. The momentum operator has the following representation in the 
position space 



So, from Eq. ( |88|1 we have 



P ^ ^ 
i dx' 



(90) 
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which has the following solution 



1 / ipx 

exp 



(91) 



V2TTh V ^ 

where the constant of integration is chosen to satisfy 

u*p{x)up{x')dp = 6{x - x'). (92) 

In the GUP scenario, to find the momentum eigenfunction in the quasiposition space, we 
write the eigenvalue equation as 

2/3 



(2/3)2/3 (^-Sifid^ + ^-dh^dl - 4//3 



1/3 



(93) 



where p is the eigenvalue of P. Now, let us take the solution in the form of Eq. f l9ip and 
fl93|) at the same time, i.e.. 



where p = f{p)- Inserting this solution in Eq. ( l93l) results in 

1 - iVS + (-2/3)1/3 (3p + ^/9p^ - 4//3 



(94) 



2/3 



(2/3)2/3 (3p + V9p2 - 4//3 



1/3 



(95) 



or 



3 



SO we have 



1 



The momentum eigenfunctions now satisfy 



exp 



W /3 3 



(96) 



(97) 



«;(O«p(0dp = / (1 - f3p') uliOupiOdp, 

sm 



2(e-n 

in agreement with Eq. ( 128|) . Finally, the plain wave ( 197|) has the maximum energy Ef 
3/(2m/3) which agrees with Eq. (H3l) . 



(98) 
(99) 
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B. Particle in a box 



As another application, let us consider a particle with mass m confined in an infinite 
one-dimensional box with length L 



V{x) 



(100) 



< X < L, 
oo elsewhere. 

The corresponding eigenf unctions should satisfy the following generalized Schrodinger equa- 
tion 



2m dx"^ 3m dx^ 18m dx^ 9m dx^ 
for < x < L and they also meet the boundary conditions '?/'n(0) = "ipniL) 



0. In 



Refs. 
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22| . the above equation is thoroughly solved to 0{f3) and its exact eigenvalues 



and eigenf unctions are found. Because of the boundary conditions, the eigenf unctions do 
not change with respect to the absence of the GUP {(3 = 0). However, the solutions exhibit 
the effect of GUP on the eigenvalues which are linear in GUP parameter j22|. These facts 
lead us to consider the following additional condition for eigenfunctions 



2m dx"^ 



< X < L, 



(102) 



where En 



2 2*2 

n n n 



2mL2 

Eq. fllOip in terms of ipn{x), i.e.. 



If this condition is also satisfied, we can write the second term in 



,x 



9x20-1) 



(x) 

^2 ) Ynyx). 



So, we have 



Hipn{x) 



En + ^PmEl + y /^^m^e^ + ^p^m^el + ■ 



i^n{x) 



Now, comparing Eqs. f llOip and fll04p shows that 



En 



En + ^PmEl + y /3'm2£3 + y/^^m^^ + ■ 

1 - + (-2)1/3 (^3^„ + v^g^^F^) 



4V3 (^3^„ + ^9^2^) 



where 7,1 = 2/3 me„. To first order of GUP parameter we have E., 



which is in agreement with the result of Ref. 



(103) 

(104) 

(105) 
(106) 

3mL^ 

22| . Moreover, because of Eq. (I102p and the 



2/3- 



1/3 



2 2fe2 

n n h 



2mL 



20 



boundary conditions, the normalized eigenf unctions are ipn{x) = -y/^sin ^ ^^^ ^ |^22]. These 
results show that in GUP* scenario there is no change in the particle in a box eigenf unctions 
but there is a positive shift in the energy spectrum which is proportional to the powers of 
/3. 

We now estimate the energy spectrum using the semiclassical scheme. For the particle in 
a box, the Wilson- Sommerfeld formula 



pdx = nh, n = l,2,---, (107) 



results in 



+ (-2/3)1/3 (sp^ + V9p2 - 4//3 

^ } -\ 1/3 

V^(2/3)2/3 3p„ + ^9pl - 4//3 



nh 

Pn = (108) 

Since the high energy momentum P depends on the low energy momentum through = 
Pn — (l/3)/3P^ dZSD, the semiclassical energy spectrum is given by 

^ 2m' 

2/3-| 2 

(109) 



It is straightforward to check that the semiclassical result f llOQp exactly coincide with the 
quantum mechanical spectrum f llOSp as well as in ordinary quantum mechanics. Therefore, 
the maximal energy of a particle in a box reads 

Emax = 7^ a- (HO) 

Imp 

C. WKB approximation 

To check the validity of the Wilson- Sommerfeld quantization rule for this modified quan- 
tum mechanics, we need to show that the zeroth-order wave function, which satisfies the 
generalized Schrodinger equation f lS^ . can be written as i/j{x) ~ exp [{i/fi) J pdx\. So let 
us take 

^(x) = e'*("), (111) 
where $(x) can be expanded as a power series in h in the semiclassical approximation, i.e., 

oo 

Hx) = -Y,^-^n{x). (112) 
n=0 
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So we have 

d'^ip{x) 



d'^ip{x) 
dx^ 



(-$'2 +2<l>")^(x), (113) 
($'4 _ 6i$'2^" _ 3$"2 _ 4^///^/ ^ (114) 



where $' indicates the derivative of $ with respect to x. Now to zeroth-order ($(a;) ~ 
$0(2;)/^) and for ^ — )■ we obtain 

+ + l^'K + =2m{E- V{x)) . (115) 

Now the comparison with Eq. flHT]) shows $q = p and consequently 



ip^x) ~ exp 



pdx 



(116) 



which is the usual zeroth-order WKB wave function obeying the Wilson-Sommerfeld quan- 
tization rule. 

VI. GENERALIZATION TO D DIMENSIONS 

We now extend the developed formalism in Sec. [IT] to D spatial dimensions. We then 
present the generalized Poisson brackets in the classical limit and study the density of 
states. 

A. Generalized Heisenberg algebra for D dimensions 

A natural generalization of the one-dimensional commutation relation ( !T6l) that preserves 
the rotational symmetry is 

where = Ylf=i PiPi- This relation implies a nonzero minimal uncertainty and a maximal 
observable momentum in each position coordinate. If the components of the momentum 
operator are assumed to be commutative 

[P,,P,]=0, (118) 
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then the Jacobi identity determines the commutation relations between the components of 
the position operator as 

X,] = '^'11 (P,X, - P,X,) , (119) 
(1 - fjP-^) 

which resuhs in a noncommutative geometric generahzation of position space. To exactly 
satisfy these commutation relations, the position and momentum operators in the momen- 
tum space representation can be written as 



P^(j){p)=p,(j){p), (120) 

ih 

Xi and Pj are now symmetric operator on the domain 5*00 with respect to the scalar product: 



X^^ip) = ^d,^4>{p)- (121) 



(^10) = / <Pp (1 - r{p)<P{p), (122) 



where = J^iLiPiPi- The identity operator is 



^^/vv. o^'^x^'- ''''' 



ipW) = V^^. (124) 



and the scalar product of momentum eigenstates is 

6^{p-p ') 
1 -/3p2 

In this representation, the components of the momentum operator are still essentially self- 
adjoint, however the components of the position operators are merely symmetric and do not 
have physical eigenstates. 

Since the commutation relations f lll7p -( fTT9l) do not break the rotational symmetry, we 
can express the generators of rotations in terms of the position and momentum operators as 

L,, = (1 - /3P') {X,P, - X,Pi) , (125) 

as the generalization of the ordinary orbital angular momentum. Now the momentum space 
representation of the generators of rotations is 

Liji^ip) = [Pidp, - Pjdp,) ip{p), (126) 

and 

[P,,L^k]=^n{5,kPJ-6,,Pk), (127) 
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[Xj, Ljk] — ih {SikXj — SijXk) , 



(128) 



[Lij, Lki] — ih (SikLji — SiiLjk + SjiLik — SjkLii) , (129) 

as well as in ordinary quantum mechanics. However, the geometry is noncommutative, 
namely 

[Xi,Xj] = (130) 

B. Density of states 

The right hand side of Eq. (fT6!) shows that the "effective" value of is P dependent. 
So the size of the unit cell in the phase space that is occupied by each quantum state can 
be also considered of as being momentum dependent. This fact changes the momentum 
dependence of the density of states and affect the calculation of the cosmological constant, 
the blackbody radiation spectrum, etc. Similar to the KMM algebra jsSi, we should first 
check that any volume of the phase space evolves such that the number of states inside it 
does not change with respect to time as the analog of the Liouville theorem. 

The Poisson brackets in classical mechanics correspond quantum mechanical commutators 

via 

hA,B]^{A,B}. (131) 

in 

Thus the classical limits of Eqs. flll7p -f lTT9|) are given by 

{^^'^.■}=(r^|7^' (132) 

{P„P,} = 0, (133) 

and the Heisenberg equations for the coordinates and momenta read (z, j run over the spatial 
dimensions and the summation convention is assumed) 

rlT-f f)M 
X, = {X„ H} = {X„ P,}— + {X„ ^.}^, (135) 

OH 

A = {P„i/} = -{X,,P,}^. (136) 
We now prove that the weighted phase space volume 

1- (3pA"d^Xd^P, (137) 
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is invariant under time evolution as the analog of the Liouville theorem. The evolution of 
Xi and Pi during an infinitesimal time interval 5t is 



X- ^ Xi + SXi, 



(138) 

(139) 



where 



SXi^ 
SP^ 



{Xi,PA^ + {x.,x,} 



dX, 



5t, 



OH 



St. 



(140) 
(141) 



After this infinitesimal evolution, the infinitesimal phase space volume is changed according 
to 



where 



d^X'd^P' 



dxi 

dX. 
dP! 



diX[, ■ ■ ■ ,X'j^,P{, ■ ■■ ,P'j- 



D) 



5. 



d{Xi., ■ ■ ■ , Xd, Pi, ■ ■ ■ , Pd) 

d6Xi dX' d5X, 



d^Xd^P. 



dSPi 



ax,' 



dX, dX, 



dPj 
dP[ 

dPi 



dP, 



d6P, 
dP, ' 



The Jacobian can be calculated to first-order in 5t as 



9{X'i, • • • , X'jj, P[, 



9{Xi, ■ ■ ■ , Xd, Pi, ■ ■ ■ , Pd) 



So we have 



.dSXj dSP, 

^ ' 9X; ^ dPi 



+ 



(142) 



(143) 



(144) 



ddXi ddPA 1 
dX. ^ 'dp J si 



d 
dXl 
d 



OH 



OH 



dX, 



{X^,P,} 



dH 



dP,'^^^''^'^dXidPj 



d 

dp 



+ 



{Xj,Pi}^ 

^ ^' dX, 



d 



■{Xi,Xj} 



dH 



+iXi,XA- — 

^ " ^^dXidXj 

d „ J dH 



_d_ 



{X,,Pi) 



dXi 

2fi{D 



{Xi,Xj} 



p. 



dx^ 

dH 



' d 

m 



2^3 



dX, 



(1-^P2) 
2pD dH 
'{l-pP^f 'dX'j' 



{X„P^} 

(1 - /3P^) 



dX, 

dXj ^ '^dPjdX. 
dH 



dxi 



- P 

2 J 



dH 
dXi 



(145) 
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which to first-order in 5t results in 



„ dH ^ 

1 ^Pi^^^t 



(146) 



Moreover 



2 



= 1-/3(P2 + 2P,5P, + ■■■), 
= l-/3(p'-2P.{X„P,}||-5t 



1 - /3P2 9Xi 
,2, , 2/3P, 



5t 



(1 - 0P' 



Therefore, to first-order in St 



D 



1 + 



1 + 



5t + - 



2(3Pi dH 
1 - pP^f dX, 



(147) 



(148) 



, (1-/3P2)' 'dX, , 
Now using Eqs. 01461) and fll48p . it is obvious that the weighted phase space volume Eq. fll37p 
is an invariant, i.e., 

(l - /3P'') d^X' d^P' = (1 - /3P2) d^X d^P. (149) 
C. The cosmological constant 

The cosmological constant can be obtained by summing over the zero-point energies of 
the harmonic oscillator's momentum states. Using the canonical form of the zero-point 
energy of each oscillator with mass m 



the sum over all momentum states per unit volume is 



(150) 



A(m) = J d^p {l - l3p^) \^-^/^~+m? 

= 2-n I dp (l — p^^Jp^ + rn?, 
Jo 



vr 



2013' 



(151) 
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where 

1 



f(x) = — VTTx(96 + 192x + 476x2 + 380x^ + 105x^) 
96 L 

-(480x2 + 720x^ + 450x^ + 105x^) cosli"^ {^/x) 1 , (152) 
and /(O) = 0. In the massless hmit we find 

A(0) = ^ = ^ [A(0)]""" , (153) 
that is ten times smaller than the massless cosmological constant predicted by the KMM 



proposal 



381]. This finite result is due to the vanishing of the density of states at high 



momenta where p = 1/ \fp plays the role of the UV cutoff. So in this scenario we do not 
need to put by hand an arbitrary scale as the UV cutoff and the cosmological constant is 
automatically rendered finite. Note that since 1/a/5 is proportional to the Planck mass 
Mpi, A(0) is too large in practice and consequently the cosmological constant problem still 
remains unsolved. However, our formulation gives the better estimation of A with respect 
to that obtained in the KMM framework. 



D. The blackbody radiation spectrum 

Because of the weight factor (1 — fSP^Y i'^ 3-dimensions, the average energy of the elec- 
tromagnetic field per unit volume at temperature T is given by 



E = — dz/ l-/3 — 



Jo \ \ c ) ) ye^^l^sT _ 

dz/M/3(z/,T), (154) 



oo 







where 



2, 3 



up{p,T)= I 1- (^) 1 «o(^^,T). (155) 



Here 



87r/iz/^ 1 

is the ordinary spectrum function and = cjh^f^. To show the effect of the minimal 
length uncertainty and the maximal momentum on the shape of the spectral function, we 
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FIG. 6. The blackbody radiation spectrum in the GUP* framework at temperature T = 0.1 T^. 



have depicted the functions 

fs(v,T) = (1 - {v/vg)^f Mv,T), (158) 
in Figs. [6] and U\ and compared them with the case of just the minimal length uncertainty 

H 

^KMM^^^y) = ^ fo{u,T), (159) 

(1 + [u/upy) 



where = c/ksy/^- As the figure shows, for small frequencies {u ^ z/^), fj3{v,T) closely 
■KMM^ However, it deviates from 



coincides with Z?'^^. However, it deviates from f™^ as the frequency increases. 



VII. CONCLUSIONS 



In this paper, we have presented a higher order generalized uncertainty principle that 
implies both a minimal length uncertainty and a maximal momentum proportional to hy/P 
and 1/a/5, respectively. The first exact representation (JTj) and (IHl) was not self-adjoint, but 
preserves the ordinary nature of the momentum operator. We found the exact eigenfunc- 
tions of the position operator and quantum mechanical and semiclassical energy spectrum of 
the harmonic oscillator in this representation and showed that the energy spectrum is also 
bounded from above. The second representation (I76|) and (177]) was formally self-adjoint and 
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resulted in the perturbative generalized Schrodinger equation in quasiposition space. We 
exactly studied the problems of the free particle and the particle in a box and showed that 
the existence of the maximal momentum Pmax = is manifest from this representation. 

We then generalized our proposal to D dimensions and found the invariant density of states. 
We showed that the blackbody radiation spectrum will be modified at higher frequencies 
and compared our results with the KMM proposal. Although the cosmological constant 
was rendered finite, the smallness of the GUP parameter resulted in a large cosmological 
constant that could not solve the cosmological constant problem. However, our calculated 
cosmological constant was a better estimation with respect to the presence of just the min- 
imal length. Here we implemented a momentum cutoff not through terms like on the 
right hand side of the commutation relations. Instead, we implemented the momentum cut- 
off through a function of P with a singularity. So the momentum space is cut into several 
sectors that decouple from each other. The sectors are separated from each other at the 
singularities of the function of P that is used. Technically, we have inequivalent irreducible 
representations of the commutation relations, one each in each sector |i41|] . This type of issue 



with the various sectors can be avoided, as it is indicated in Ref. j42|. Applied to our case, 
the trick would be to write the right hand side of the commutation relation not as a fraction 
but instead to expand it out as a geometric series. It has a finite radius of convergence and 
that rules out all representations beyond the singularity. 
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Appendix: Harmonic oscillator spectrum in the KMM framework 

In the context of the KMM proposal, the total energy in terms of ordinary variables is 
given by 

E= I^^ + -muj^l + fip'? x\ (A.l) 
Now the Wilson- Sommerfeld integral can be written as 

xdp = — r "{^^'f dp = (n+l) h, (A.2) 



muj 1 + \ 2 

where z = \/2mE. So the semiclassical energy spectrum is given by 

^) = -linoo + nuj(^n + 1^ (l + I) + llhoon^ (A.3) 

where 7 = (imhuj. This result agrees (up to a constant) with the exact solution to first order 
of the GUP parameter [7] 

ET"' = noo{n + ( + + 7/2) + linun' 

= i?f^) + ^7^a; + 0(7^), (A.4) 

and gives the correct dependence behavior. An alternative derivation of Eq. (1A.3P is also 



presented in Ref. [251]. 



ACKNOWLEDGMENTS 



I am very grateful to Achim Kempf for fruitful discussions and suggestions. 



[1] G. Veneziano, Europhys. Lett. 2, 199 (1986). 
[2] E. Witten, Phys. Today 49, 24 (1996). 

[3] D. Amati, M. Ciafaloni, G. Veneziano, Phys. Lett. B 216, 41 (1989); Nucl. Phys. B 347, 550 

(1990); Nucl. Phys. B 403, 707 (1993). 
[4] K. Konishi, G. Paffuti, P. Provero, Phys. Lett. B 234, 276 (1990). 
[5] L.J. Garay Int. J. Mod. Phys. A 10, 145 (1995). 
[6] M. Maggiore, Phys. Lett. B 319, 83 (1993). 

30 



[7] A. Kempf, G. Mangano, R.B. Mann, Phys. Rev. D 52, 1108 (1995). 

[8] A. Kempf, G. Mangano, Phys. Rev. D 55, 7909 (1997). 

[9] M. Maggiore, Phys. Lett. B 304, 65 (1993). 

[10] F. ScardigU, Phys. Lett. B 452, 39 (1999). 



[11] M. Maggiore, Phys. Rev. D 49, 5182 (1994), (a?Xiv:hep-th/93d5T63 



[12] S. Hossenfelder et al, Phys. Lett. B 575, 85 (2003), arXiv:hep-th/0305262 



[13] C. Bambi, F.R. Urban, Class. Quantum Grav. 25, 095006 f2008). [a?Xrv:0709.1965l 

[14] K. Nozari and B. Fazlpour, Gen. Relativ. Grav. 38, 1661 (2006). 

[15] R. Banerjee and S. Ghosh, Phys. Lett. B 688, 224 (2010). la?Xiv:1002. 23021 

[16] P. Pedram, Int. J. Mod. Phys. D 19, 2003 (2010). 

[17] K. Nozari and P. Pedram, Europhys. Lett. 92, 50013 f2010). [arXiv:1011. 56731 
[18] B. Vakih, Phys. Rev. D 77, 044023 (2008). 

[19] M.V. Battisti and G. Montani, Phys. Rev. D 77, 023518 (2008). 
[20] B. Vakih and H.R. Sepangi, Phys. Lett. B 651, 79 (2007). 
[21] M.V. Battisti and G. Montani, Phys. Lett. B 656, 96 (2007). 
[22] K. Nozari, T. Azizi, Gen. Relativ. Gravit. 38, 735 (2006). 

[23] S. Das and E.G. Vagenas, Phys. Rev. Lett. 101, 221301 f2008). la?Xrv:0810.5333[ 

[24] K. Nouicer, Phys. Lett. B 646, 63 (2007). 

[25] P. Pedram, Phys. Rev. D 85, 024016 f2012). [a?Xrv:1112.23271 

[26] P. Pedram, Phys. Lett. B (2012), DOL10.1016/j.physletb. 2012.03.015. 



[27] J. Magueijo and L. Smolin, Phys. Rev. Lett. 88, 190403 (2002), |arXiv:hep-th/0112090 



[28] J. Magueijo and L. Smolin, Phys. Rev. D 71, 026010 (2005), |arXiv:hep-th/0401087; 

[29] J.L. Gortes and J. Gamboa, Phys. Rev. D 71, 065015 (2005), [a^v:hep-th/0405285 

[30] P. Pedram, Europhys. Lett. 89, 50008 (2010). la?Xrv:1003.2769l 

[31] P. Pedram, K. Nozari, and S.H. Taheri, JHEP 03, 093 (2011). 

[32] A.F. Ah, S. Das, and E.G. Vagenas, Phys. Lett. B 678, 497 (2009). 

[33] S. Das, E.G. Vagenas, and A.F. Ah, Phys. Lett. B 690, 407 (2010). 

[34] A.F. Ah, S. Das, and E.G. Vagenas, Phys. Rev. D 84, 044013 (2011). 



[35] A. Kempf, arXiv:hep-th/9603115 



[36] P. Pedram, M. Mirzaei, and S.S. Gousheh, Mol. Phys. 108, 1949 (2010). 



31 



[37] We have used the relation Sh ^^ = 2ri'(^l- ^yl-2rf-^{n + l/2)j . 

[38] L.N. Chang, D. Minic, N. Okamura, T. Takeuchi, Phys. Rev. D 65, 125028 (2002). 

[39] N.I. Akhiczer and I.M. Glazman, Theory of Linear Operators in Hilbert Space (Dover, New 
York, 1993). 

[40] G. Bonneau, J. Faraut, and G. Valent, Am. J. Phys. 69, 322 (2001). 

[41] Private communication with Achim Kempf. 

[42] A. Kempf, Phys. Rev. Lett. 92 (2004) 221301. 



32 



